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SUMMARY 


Recent advances in electronics technology along with the advent of low cost 
multi-channel Fast Fourier analyzers have now made it practical to use higher 
order central difference formulae to measure power flow in one dimensional and 
two dimensional structures . The method discussed in this paper uses five point 
differencing for the spatial derivatives in one dimension and a thirteen point 
difference pattern for the spat lal derivatives in two dimensional plates and 
shells . It is assumed that the measuring transducers are accelerometers . 

An analytical study of the higher order differencing method and the conven- 
t ional two accelerometer method was performed here as a preliminary to the appli- 
cation of these methods to actual aircraft structures . Some classical problems 
were analyzed in order to simulate and compare the performance of the two methods 
under near field measurement condition? . Near field conditions analyzed in this 
study include examples of power flows near simple sources and simple boundaries 
The estimates produced by the two methods were compared to the exact solut ion in 
each example . This paper presents the theory and selected results of the study . 
The results indicate that the bias errors of the two accelerometer method under 
near f leld measurement condit ions may be much larger than previous studies have 
suggested 

INTRODUCTION 

The successful realization of a diagnostic measurement device for measux ing 
power flow would greatly aid design engineers in their efforts to reduce the air- 
craf t interior noise and vibration of both fixed wing and rotary wing aircraft . 
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For example , accurst* power flow measurements would allow designers to deter- 
mine the critical paths of the vibrational energy transmission through the var- 
ious structural members, verify computer predictions produced by sophisticated 
analytical models , and determine the relative level of effectiveness of various 

treatment methods used in reducing the vibrational power transmission in a struc- 
ture . 

In order for a power flow measurement method to be utilized to its fullest 

potential on aircraft structures it must meet the following operational require- 
ments : 

( 1 ) The method must be accurate when measuring in evanescent or reverberant flow 
fields . 

(2) The method must be applicable to a variety of aircraft structures with 
different material properties and aspect ratios. 

(3) The method must be applicable over a wide range of frequencies . 

Nearly all of the research performed in the past fifteen years on experimen- 
tal measurement methods for measuring power flow in structures has focused on the 
two accelerometer method. 1 '■ In order to circumvent the difficulties of keeping 
track of the multiplicity of terms required to measure power flow, the two ac- 
celerometer method utilizes several simplifying assumptions, sometimes called 
Noxseux ‘a assumptions' , which are based on the premise that the shear wave compo- 
nent and the bending wave component are equal under "free field" measurement con- 
ditions (A "free field" contains no power sources or boundaries in the vicinity 
of the transducers ) 

As a result of the use of these assumptions, it was known at the outset that 
the data obtained by Noiseux * s method would be inaccurate near power sources 
or boundaries and that only an estimate of the total power flow with no shear- 
bending- twist component breakdown would be possible. The method has remained 
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popular, however, ainca it requires only two transducers for measurements in a 
/ ] 
one dimenional structure (four transducers for a two dimensional structure) and 

1 

only a dual channel Fast Fourier analyzer to perform the spectral analysis of 
the data. Largely due to the simplicity of the measurement system, this method 
appears to have become widely accepted for measuring power flow in simp] e beams 
and plates . The two principal difficulties with the method have been identified 
by researchers as inaccuracies while measuring under near-field conditions 4 ’ 6 , 

i 

and signal-to-noise problems associated with highly reveberant measurement 
conditions . 45 

Recent advances in electronics technology along with the advent of low cost 
multi-channel Fast Fourier analyzers have now made it possible to use mere so- 
phist icated methods to measure power flow in one dimensional and two dimensional 
structures Extremely low weight (0.3 gram) miniature piezoelectric accelerom - 

t 

e t er 3 can now be used in conjuctionwith a multi-channel Fa3t F our i er analyzer i 

to perform power f low measurement s . The method proposed in this paper (based in 
part on Pavic’s earlier work-) performs „ direct finite difference approxima- 
tion of the spatial derivatives by ut i 1 izing u 13 accelerometer computational 
molecule for the two dimensional problem ( 5 accelerometers for the one dimen- 
sional problem) The primary advantage of this method i s that its results remain 
accurate near power sources and boundaries and it provides a shear-bending-twist 
component breakdown of the total power flow (With the large number of power 
sources and the hundreds of structural junctures and boundaries in an aircraft , 
these improvements are believed to be essential ) 

An a first step in applying this new method to aircraft structures , an 
analyt ical study of the new method and the older two accelerometer method was 
performed in order to simulate and compare the performance of the two methods 
under near field measurement condit ions . Measurement conditions analyzed in 



thi» study include simulated measurements near 1 


ine forces, line moments, point 


P'’ i " e Iniaral quadrupoles , and boundari.a »hich p».,.„ 

varioua to«bi...,o„, af «... , translational .tiffn... . «... oi ia „ n ^ 

and rotational stiff.... prop.rtia. Th. pr.drct.d r.sult. of th. t.o „. th „d. 
-er. compared to the exact solut ion m each case . This _ 


ln each • Thrs paper outlines the theory 

dtiliz.d in thi. study and pr.aanta th. coitput.d r.ault, for th. .inul.t.dn.ar 
field measurement conditions 


SYMBOLS 


amplitude vector of a right traveling wave at position , |) 


amplitude v.ctoi of a left traveling wave at position .r () 


J ) l »^or of a right traveling wave at position r 


amplitude vector of a left traveling wave at position r 
b ending stiffness of the plate given by the expression Fh\'[ I 2 ( 1 ,,2)j 

’ < : arbitrary amplitude coefficients used in the general solution of the 


governing equation 


reflection matrix coefficient used to represe.it the translational 


stiffness and inertia 


r.f lsct ion matrix co.ffici.nt un.d to r.pr.s.ut th. rotational 


and inert ia 


base of the natural logarithms 


Young’ s modulus (modulus of elasticity) 


f requericy 


propagation matrix for a right traveling wave 


propagation matrix for a left traveling wave 


amplitude (source strength) of a point force input 


auto spectrum of acceleration signal i , 2 , 3 . 13) 


cross spectrum of acceleration signals <• and . f (<»,.? = j 2,3 13) 


© 
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I> plat# thickness 

Hl"{ ) Oth order Hankel function of the first kind 
Hl ~ >( ) Oth order Hankel function of the second kind 
H\'\ ) 1st order Hankel function of the second kind 
) 2nd order Hankel function of the second kind 
/ area moment of inertia about the neutral axis 

j/| identity matrix 

imaginary part of the quantity inside the braces 
/ square root of -1 

A vavenumber 

A„ translational stiffness of the boundary 

i„ rotational stiffness of the boundary 

l\ dimensionless translational stiffness 

l\„ dimensionless rotational stiffness 

m mass per unit width for a beam , mass per unit area for a plate . 

m n mass per unit length of the boundary 

n , rotational inertia per unit length of the boundary 

\l h bending moment per unit length in the direction of wave propagation 

1/,,. bending moment per unit length perpendicular to the direction of 

wave propagation 

.1/ amplitude ( source strength) of a point moment input 

Mi twi st ing moment per unit length in the plate 

P loadi ng f unct ion 

/' amplitude of the loading f unct ion 

( ) shear per unit length in the plate 

Q amplitude ( source 5 Irength) of a quadrapole input 

, radial distance from the source 
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ref lection matrix 



real part of the quantity inside the braces 


cartesian coordinate 
cartesian coordinate 

real power per unit length flowing in the x direction 
real power per unit length flowing in the y direction 
bending component of the real intensity in x direction 
shear component of the real intensity in x direction 
twist component of the real intensity m x direction 
bending component of the real intensity in y direction 
shear component of the real intensity in y direction 
twist component of the real intensity in y direction 
denotes partial differentiation 
Dirac delta function 

transducer (finite difference) spacing 

transverse d i splacement of the beam or plate 

transverse displacement at location i on the on the beam 

or piste (1 = 1, 2, B , •• 13) 

solution to the modified Bessel ' s equation 

solution to Bessel's equat ion 

slope in the direction of wave propagation 

angle of twist 

Poisson ’ s rat 10 

translational damping ratio for the boundary 
rotational damping ratio for the boundary 
ratio of ci rcumf nrence to the diameter of a circle 
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time averaged complex intensity in ^ha x direction 


azimuthal angle 


radian driving frequency 


Laplacian operate 


Biharmonic operator 


denotes a time average 


approximately equal to 


differentiation with respect to time 


twice differentiation with respect to time 


REVIEW OF THE FREE FIELD METHOD 


The time averaged intensity (power flow per unit length) of transverse waves 
flowing in the .r direction through a two dimensional structure is given by the 


equat ion 
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where the individual terms in equation ( 1 ) are given by 
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Noiseux 1 has shown that equation (1) can be simplified considerably when the 
power flow takes place under free field conditions . In reference 1 , it is shown 
that if there are no power sources or discontinuities in the properties of the 
media (boundaries ) nearby , then the following relationship holds: 


M h ft h | M,», -f 


A Mb f Mb, 


Further, it is shown in reference 1 that the approximation given by equatic 
(8) can be easily obtained from measurable quantities under free field measure- 
ments conditions from the following relation: 


Mb + Mb-, I. 

<T /«)" : 


//V-r/ M-’,,. 


Substituting the approximation of equation (9) into equation (1) the re- 
sulting equation for the total power flow is 

Nj. ■ Qij ■ , < Ilk 1 ■ >)ft h -i ( |f 

Since the total power flow is shared equally by the shear and bending com- 
ponents under free field conditions , equation (10) can be reduced further to the 
following simple result : 

MU- nil. . m 


Equation (11) has been useu by a number of researchers to obt 


am an estimate 


of the power flowing in simple beam and plate structures' " . The most populi 
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experimental implementation of equation (11) utilizes two accelerometers and the 
following finite difference approximations : 


>1 - I'/l i 


( 12 ) 


= \vi m}/&- 

Substituting equations (12) and (13) into equation (11) the result is 
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Expressing equation (14) in terms of acceleration rather than displacement and 
velocity , converting to the frequency domain and taking the real part of the 
result , the equation for intensity is 


in.' 2 in* , ,,, , 

• /?n{kV|| (V>j t ( *2 1 0*12]} — ^ j * *** l^ r 2l i ■ 


and since A - " 1 u/-m /I the result can be written as 
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Thus , equation ( 16) can be used as an estimate of the power flowing in a beam 
or plate using the two-accelerometer method 

THE DIRECT FINITE DIFFERENCE APPROACH 

An alternate approach to the problem of neasuring power flow in structures 
is to perform direct finite difference approximations of equations (2) through 
(7) by using a more sophisticated computational molecule Consider the compu- 
tational molecule shown in figure 1 This computational molecule consists of 
an array of 13 measurement transducers (accelerometers) arranged in a symmetric 
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fashion on a two dimensional surface (e . g . a plate) . If each transducer in the 
molecule consists of a 0 . 3 gram accelerometer, the combined mass load on the sur- 
face of the structure amounts to 3 . 9 grams . (The mass load per unit area will de- 
pend on the tranducer spacing . ) Other computational molecules with fewer trans- 
ducers (e.g. 10 transducers 2 ) can be used in the finite difference approach . The 

13 transducer molecule, however, was chosen for study since it provides estimates 
of the power flow in the r and y directions at the same location in a 2 dimensional 
structure . 

Utilizing the computational molecule of figure 1 , the following central 
difference approximations can be made : 
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These central difference approximations may be rewritten in terms of ac- 
celerations instead of displacements and substituted into the appropriate terms 
into equation ( 1 ) Converting this form of equation (1 ) to the frequency domain 
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and taking the real part of the result , it can be shown that the shear , bending , < 

and twist components of the power flow are given respectively by 
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Recalling the 
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following properties of cross spectra 
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0 , (2(5) 

(27) 

it 1 a seen that equations (23) through (25) require a total of 21 independent 
cross channel measurement pairs in order to compute estimates of the real part of 
the shear . bending, and twist components of the power flow . Thus , the experimen- 
tal imp 1 ementat ion of the computational molecule of figure 1 requires at least 2 
passes with most 16 channel FFT analyzers . or 3 passes with most 8 channel FFT an- 
alyzers . 

Finite difference approximat ions for the sheer , bending , and twist compo- 
nents of the power flow in the 7 direction can be obtained in a similar fashion. 

The results of the derivations are given by the equations 
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Equations (28) through (30) indicate that 13 additional independent cross 
channel measurement pairs are required to obtain the component s of the real power 
flow in the y direction. Thus, the implementation of the computational molecule 
of figure 1 requires a total of 34 independent cross spectral measurements to 
obtain the six components of power flow in e r and y directions . In practice, 
this would require 3 passes with most 12 or 16 channel FFT analyzers, or 5 passes 
with most 8 channel analyzers. 

It should also be noted that if the geometrical characteristics of the 
structure under investigation are 1 dimensional , the direct finite difference 
approach simplifys considerably For example , if the structure is a simple beam, 
the equation for the power flow is given by 

ll r Q’l -r ' - Mh^b •< » 

For the case of the 1 dimensional structure, the central finite difference 
approximations for >) and for remain the same (see equations (17) and (18)) . The 
formulae for the shear and bending moment simplify to the expressions 


. <^7 , i 'h 2 i/i 1 27m 7 n I 
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Thus, it is seen that only 5 of the 13 transducers are required in this case. 
Further analysis shows that the equations for the real parts of the shear and 
bending components of the power flow are given by 
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Equations (34) and (35) indicate that only 5 independent cross spectral 
measurements must be performed in order to obtain estimrtes of the real part of 
the shear and bending components of the power flow in a beam. This requires a 
single pass with most multichannel FFT snalyr.ors . 


SIMULATED MEASUREMENTS 

A series of analytical simulations of the two accelerometer and the direct 
finite difference measurement methods were performed in order to quantify and 
compare the accuracy of the two methods under near field measurement conditions . 
In the subsections that follow, selected results of the computer simulations 
of the two methods are presented for the cases of measurements in a plate near 
to a line force, a line moment , a point force, apoint moment , and a quadrupole . 
Secondly , selected results are presented for simulated measurements near to 
boundaries which possess various combinations of mass , translational stiffness , 
rot at l onal inertia, and rotational stiff ness proper t ies . (The geometrical 
arrangement of the i nput fore ing f unc t ion , the transducers , and the boundary 
are shown in f igure 2 .) A 032 cmthickAA2024 aluminum plate was chosen as the 
test vehicle for the simulations since it is representative of the materials and 
structural aspect ratios used in aircraft construction. The material properties 
of the plate are presented in Table I 
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Each data figure presented in the subsections that follow consists of parts 
(a) and (b) . Part (a) shows the relative error of the two measurement methods in 
dB (relative to the known exact solution) over the 0 - 1000 Hz frequency range. 

Part (b) shows the intensity of the power flow in dB (relative to 1 Watt/meter) 
as estimated by the two measurement methods for a forcing function which ha9 a 
uniform input over the 0 - 1000 Hz frequency range . In each case , the two figures 
viewed together provide some insights into the significance of the error The 
results produeedby the two accelerometer method are denoted by the words FREE 
FIELD and the results produced by the direct finite dif ference approach are 
denoted by the word DIRECT in both the figure legends and in the discussions that 
follow . The 0 - 1000 Hz frequency range was chosen for the analysis since it is of 
primary importance in propeller driven aircraft 

The exact solution for the structural intensity flowing in the plate for 
each flow field under examination was obtained in the following manner: The ex- 
act solutions for the transverse displacement , transverse velocity , angular vo- 

locity . rate of twist . shear , bending moment , and twisting moment as a function of 
posit ion on the plate were derived (Summaries of these analytical derivations 
are presented in Appendices I through III ) The results of these derivations were 
incorporated in the computer codes and were used in conjunction with equation ( 1 ) 
to calculate the exact amounts of structural intensity flowing in shear , bending, 
and twist at any location in the plate 

The two measurement methods were simulated fi r each flow field under inves- 
tigation as follows The analytical result for the transverse displacement as a 
function of position on plate was utilized to compute the acceleration that would 
be expenen ed by each individual accelerometer m the transducer array . (The 

computed transverse displacement was multiplied by m the frequency domair 

to obtain the acceleration in each case ) The computed accelerat ions were then 
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Each data figure presented in the subsections that follow consists of par ts 
(a) and (b) . Part (a) shows the relative error of the two measurement methods in 
dB (relative to the known exact solution) over tho 0 - 1000 Hz frequency range 
Part (b) shows the intensity of the power flow m dB (relative to 1 Watt/meter) 
as estimated by the two measurement methods for a forcing function which has a 
'mm uniform input over the 0 - 1000 Hz frequency range In each case , the two figures 

^ viewed together provide some insights into the significance of the error . The 

0* results produced by the two accelerometer method are denoted by the words FREE 

FIELD and the results produced by the direct f inite difference approach are 
denoted by the word DIRECT in both the figure legends and in the discussions that 
follow . The 0 - 1000 Hz frequency range was chosen for the analysis since it is of 
primary importance in propeller driven aircraft 

The exact solution for the structural intensity flowing in the plate for 



substituted into the finite difference approximations utilized by each measure- 
ment method. Thus . a 4 transducer computational molecule was used in conjunction 
with equations (12) and (13) to simulate the two accelerometer method while a 13 
transducer computational molecule was used in conjunction with equations (18) 
through (23) to simulate the direct finite difference approach. 

A 2 cm transducer spacing was selected for the two accelerometer method 
while a 1 cm transducer spacing was selected for the direct finite difference ap- 
proach. It is acknowledged at the outset that these choices provide the direct 
approach with a slight advantage over the two accelerometer approach in the high 
frequency ranges (beyond the scope of analysis of this paper) where finite dif- 
ference error is important In the low frequency ranges where near field error 
is important, however, the larger transducer spacing for the two accelerometer 

probe should, if anything, provide it with a slight advantage over the direct ap- 
proach 

measurements near a line force 

Figures 3 and 4 show typical results obtained for simulated measurements 
m a plate near to a line force input A summary of the analysis used to generate 

these results .s given in Appendix I The amplitude of the line force disturbance 
was fixed at 10.0 N/m 

For the case shown in figure 3. the geometric center of each measurement 
probe is located at a distance of 10 cm from the line force discontinuity (31.5 
times the thickness of the plate) Figure 3(a) indicates that the near field er- 
ror for the FREE FIELD method is less than 1 dB everywhere except in the frequency 
range below 50 Hz The near field error of the DIRECT method is seen to be negli- 
gible in ail frequency ranges Fig Ure 3(b) adds additional insight into figure 
3(a) since it indicates that the frequency range m which the FREE FIELD method 
experiences the largest errors coincides with the frequency range of maximum 
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poser f los for a uniform input force . Thus , the error in the estimates produced 
by the FREE FIELD method can be significant in some cases (Note : References 4 , 

5 . and 6 present results similar to those shown in figure 3(a) for the case of a 
line force input , but do not include the accompanying plot of the intensity of the 
power flow ) 

Figure 4 shows the results obtained when the measurement probes are located 

6 cm from the line force discontinuity (15.75 times the plate thickness) . These 
results indicate that the magnitude of the error in the intensity estimates 
produced by the FREE FIELD method increases in the lower frequency ranges . The 
est imates produced by the DIRECT method remains very accurate , however . Also , 

f rom comparison of figures 3 and 4 , it can be seen that the errors produced by 
the free field approximations will change sign in some frequency ranges when 

the probe i s moved closer to the input . Thus , the FREE FIELD method may either i 

i 

overestimate or underestimate the true intensity depending on the location of the 
measurement probe . 

MEASUREMENTS NEAR A LINE MOMENT 

F igure 5 shows typical results obtained for s imulat ed measurement s in a 
plate near to a 1 me moment input . The analysis usad to generate these results is 
summarized in Appendix I . The amplitude of the line moment disturbance was fixed 
at 0.20 Nm/m . For the cass shown in figure 5, each measurement probe is located 10 
cm from the line moment discontinuity . In this case, the near field error in the 
eat imates produced by the FREE FIELD method is seen to be con si stent ly negative 
(the true intensity is underestimated) and the errors are larger in magnitude 
than those experienced when measuring near a line force (see figure 4) . The 
near field error of the est imat es produced by the DIRECT method in f igure 5 is 
negl i g lb] e by compar i son . These results are in sharp contrast to the results of 
f igures 3 and 4 where the error curve produced by the FREE FIELD method "hovers" 
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about the axis of 0 error. In fairness to the FREE FIELD approach, however, 
it should be noted that the frequency range in which large error occurs in its 
estimates coincides with the frequency range of minimum power flow for a uniform 
line moment input . 

MEASUREMENTS NEAR A POINT FORCE 

The two methods of power flow measurement were also simulated for the case of 
measurement near a uniform point force input to an infinite plate . The magnitude 
of the point driving force input was fixed at 10 . 0 N . Qualitatively , the results 
obtained for these simulations were very similar to those obtained for a line 
force input (See figures 3 and 4 . ) . This is to be expected since the solutions to 
this type of flow field (summarized in Appendix II) must obey the radial symmetry 
of the problem. Therefore, the results obtained should exhibit similar one 
dimensional type behavior. 

The simulations for the point force did produce one notable addition to 
the previous body of Know lege , however r or the case shown in figure 6 , each 
measurement probe is located 5 cm from the point force discontinuity . In figure 
6(a) i the error curve for the FREE FIELD method does not cross the ax i s of 0 
® rror • but underestimates the true intensity over the entire frequency range of 
interest (In other cases where the measurement probes are 10 cm or further from 
the source, the error curve of the FREE FIELD method crosses the 0 error axis. ) 

Thus , unlike the results obtained for the line force, the error curve of the FREE 
FIELD method does not necessarily "hover' about the axis of 0 error . Since the 
amount of power input to the plate is uniform over the entire frequency range in 
this case (see figure 6(b)) , the large near field error introduced by the FREE 
■■ IlLC method at low frequency could significantly degrade the accuracy of the 
estimate of the total power flow Also note that the noar field error experienced 
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by th« DIRECT method of measurement is once again insignif icant in comparison to 
that experienced by the FREE FIELD method . 
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MEASUREMENTS NEAR A POINT MOMENT 

Figures 7 and 8 show typical result s obtained (using the analysis summarized 
in Appendix II ) for simulated measurement sinaplate near to a point moment 
input . The amplitude of the point moment was f ixed at 0 . 20 Nm . Unlike the result s 
presented for the previous cases, the flow field in this case possesses a strong 
directional character. 

For example figure 7 shows the result s of the simulations when each mea- 
surement probe is located 10 cm from the point moment and aligned along the dipole 
axis ('/> 0 degrees) (The flow field of a point moment can be mathematically 

represented as a force couple or dipole.) In this case, the near field error in 
the estimates produced by the FREE FIELD method is consistently negat ive , but l s 
unacceptably large only at frequencies below 200 Hz . The DIRECT method is seen to 
produce very accurate estimates by comparison 

As the measurement probes are moved either radially inward (closer to the 
source ) , or moved circumferentially in the azimuthal angle •:> , the performance 
of the FREE FIELD method is seriously degraded . For example , f igure 8 shows the 
results of the simulat ions when each measurement probe is located at an azimuthal 
angle O of 80 degrees The measurement probes are again directed at the point 
moment source and locatedat adistanceof 10 cm away . Thus, tke results of figure 
8 are for the measurement of intensity along an axis that is almost in pure twist. 
These results indicate that the estimates produced by the FREE FIELD method along 
this axis are so inaccurate that no useful information can be obtained from them 
In contrast, the estimates produced by the DIRECT method remain accurate across 
the ent i re 0 - 1000 Hz frequency range 
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measurements near a quadrupole source 

Figures 9 and 10 show typical results obtained (using the analysis summa- 
rized in Appendix II) for simulated measurements in a plate near to a lateral 
quadrupole input The amplitude of the quadrupole was fixed at 0.02 Hm- . The 
case of the lateral quadrupole is of interest since its flow field closely resen. 
bias that of a point force input located near the corner of a semi- inf imte plate 
which is simply supported along its two edges 

Figure 9 shows the results of the simulations when each of the measurement 
probes is directed at the quadrupole source lccated 10 cm away . In this case, the 
azimuthal angle >p was fixed at 45 degrees This should be the best case sceneno 
for FREE FIELD measurement accuracy since the axis along the 45 degree azimuth 
is the axis of minimum twist . The near field error in th« estimates produced by 
the FREE FIELD method is consistently negative and significant at frequencies 
below 300 Hz . The DIRECT method is seen to produce very accurate estimates over 
the entire frequency range 

Qnce again, as the measurement probes are moved either radially inward, or 
moved circumf erent lally , the performance of the FREE FIELD method is seriously 
degraded . Figure 10 shows the results of the simulat ions when each of the mea- 
surement probes is located at an azimuthal angle O of 10 degrees The probes are 
again directed at the quadrupole source and located at a distance of 10 cm away - 
The estimates produced hy the FREE FIELD method along this axis are grossly inac- 
curate The estimates produced by the DIRECT method, however , remain accurate 
across the entire frequency range 

MEASUREMENTS NEAR A LINE BOUNDARY 
WITH MASS AND TRANSLATIONAL STIFFNESS 
Figures 1 1 and 12 show typical results obtained from the simulations for 
measurements in a plate near to a line boundary which possesses a translational 
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stiffness A-„ of 1.0 * I0 T H/m/m, amass »»,, of i.Okg/m, and a viscous critical damp- 
i-- ratio („ of 0 0 I . A summary of the analysis used to generate the results is 
given in Appendix III. The analysis utilizes a uniform line force input of 10.0 
N/m amplitude located 1 m away from the boundary to generate incident travel- 
mg waves on the boundary (see figure 2) The numerical values chosen to char- 
acterize the boundary and the input forcing function are somewhat arbitrary but 
representative of typical values that might be encountered in existing aircraft 
type structures . The large distance between the source and the boundary insures 
that the boundary does not experience any near field effects associated with the 
source . Thus, the boundary ’’sees " only incident propagating waves . 

For the case shown in figure 11 , the measurement probes are located 10 cm 
away from the boundary Figure 11(a) indicates that the near field error for 
the FREE FIELD method is less than 1 dB everywhere except in the frequency range 
helow 300 Hz. The DIRECT method is seen to be accurate over the entire frequency 
range Figure 11(b) indicates that the frequency ranges of largest error and 
maximum power flow coincide for the FREE FIELD method Thus, the error m the 
est mates produced by the FREE FIELD method can be significant . Also note that 
very little power flows through the boundary in the vicinity of 300 Hz. This 
phenomenon is caused by the strong impedance mismatch effects of the boundary 
at this frequency It can be shown that the exact location (in frequency) of the 
impedance mismatch is dependent primarily on the value of the stiffness of the 
boundary A ,,. K 

Figure 12 shows the results obtained when t ie measurement probes are located 
5 cm from the boundary Figure 12(a) indicates that the magmtudeof the error in 
the intensity estimates produced by the FREE FIELD method increases as the trans- 
ducers are moved closer to the boundary The estimates produced by the DIRECT 
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method remain very accurate , however. Figure 12(b) is plot of the absolute val- 
ues of the intensity estimates for the two methods . It is important to note that 
the FREE FIELD method produces negat ive intensity values in the 175 - 275 Hz fre- 
quency range . Thus , the FREE FIELD method indicates that the power is flowing out 
of the boundary rather than through the boundary in thi s frequency range . 
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MEASUREMENTS NEAR A LINE BOUNDARY WITH 
ROTATIONAL STIFFNESS AND INERTIA 

Figures 13 and 14 show typical resulcs obtained (using the analysis of 
Appendix III) from the simulations of measurements in a plate near to a line 
boundary which possesses a rotational stiffness A w of 1 .0 * 10 1 N-m/m, a rotational 
mass moment of inert la in,/ of fUHH kg-m-'/m, and a viscous critical damping ratio 
ofOOi. The analysis utilizes a uniform line force input of 10.0 N/m amplitude 
located 1 m away from the boundary to generate incident traveling waves on the 
boundary (see figure 2) Once again, the numerical values chosen to characterize 
the boundary and the input forcing function are arbitrary but realistic values . 
The large distance between the source and the boundary insures that the boundary 
does not experience any near field effects associated with the source. 

For the case shown inf igure 13, the measurement probes are located 10 cm 
away from the boundary . Figure 13(a) indicates that the near field error for 
the FREE FIcLD method is acceptable over the enti re frequency range. The DIRECT 
method is seen to be slightly more accurate than the FREE FIELD method. Also 
note that very little power flows through the boundary inthe vicinity of 850 Hz . 
Thus, the impedance mi smatch phenomenon occurs at a much higher frequency in this 
case It can be shown that the location o* the impedance mismatch is dependent 
primarily on the value of the rotational inertia of the boundary 

Figure 14 shows the results obtained when the measurement probes are located 
5 cm from the boundary Figure 14(a) indicates that the magnitude of the erroi 
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in the intensity estimates produced by the FREE FIELD method greatly increases 
at frequencies above 500 Hz. Tho estimates produced by the DIRECT method remain 
very accurate, however. Figure 14(b) is plot of the absolute values of the in- 
tensity estimates for the two methods . Mote that the FREE FIELD method produces 
negative intensity values in the vicinity of 800 Hz indicating that the power is 
flowing out of rather than through the boundary in this range Therefore, the ac- 
curacy of the FREE FIELD method can not necessarily be guaranteed in the higher 
frequency ranges , but may depend on the location of the transducers and the char- 
acteristics of the boundary 

MEASUREMENTS NEAR A LINE BOUNDARY WITH 
TRANSLATIONAL AND ROTATIONAL PROPERTIES 
Figures 15 and 16 show typical results obtained (using the analysis of Ap- 
pendix III) for the simulations of measurements in a plate near to a line bound- 
ary which possesses both the translational and rotational stiffness and inertia 
properties discussed in the two previous subsections 

( In ► MR N/m/m, »i,. 10 kg/m. 0.01 , 

l.o. Ill 1 N-m/rn, n>„ 0.001 kg-m-/m,and .01 ). 

The analysis utilizes a uniform line force input of 10.0 N/m amplitude located 1 m 
away from the boundary to generate incident traveling waves on the boundary (see 

figure 7 . ) 

Figure 15 shows the results obtained when the measurement probes are located 
iO cm away from the boundary These results bear a strong resemblance to those 
obtained in figure 1 1 . Figure 15(a) indicates that the near field error for the 
FREE FIELD method is less than 1 dB everywhere except in the frequency range of 
maximum power flow. The DIRECT method is seen to be accurate over the entire 
frequency range Figure 15(b) suggests that the added rotational propert les of 
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the boundary cause a secondary impedance mismatch at some frequency above 1000 
Hr. 

Figure 16 shows the results obtained when the measurement probes are located 
5 cm from the boundary The error curve of figure 16(a) indicates that the inten- 
sity estimates produced by the FREE FIELD method are inaccurate over the entire 0- 
1000 Hz frequency range . The estimates produced by the DIRECT method remain ac- 
curate , however . Figure 16(b) shows that the intensity estimate t> produced by the 

/ 

FREE FIELD method have the wrong sign near the first impedance mismatch (200-300 
Hz) and deteriorate rapidly as they approach the second impedance mismatch (above 
1000 Hz) . Thus , it may be cone luded that the magnitude and character of the n< ar 
field error produced by the FREE FIELD method may depend on the particular combi - 
nation of the translational and rotational properties of a boundary in the vicin- 
ity of the measurement probe 

MEASUREMENTS UNDER SEVERE, 

COMBINED NEAR FIELD CONDITIONS 
Figures 17 shows typical results obtained (using the analysis of Appendix 
III) for the simulations of measurements in a plate with the measurement probes 
"sandwiched" between a line force and a line boundary. The boundary in this case 
possesses the proper t les discussed previously 

( k\ t I II » III 7 N/m/m, u>., 1.0 kg/m, 11.0! , 

1.0 • 10* N-m/m, >>>,> 0.00| kg-m J /m, and .01 ) , 

and is located 5 cm from the measurement probes. A line force input of uniform 
10.0 N/m magnitude is posit loned lOcmfromthe boundary Thus , the distance 
between the input and the geometric center of the transducers is also 5 cm. 
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Under these conditions , the transducers are affected by the propagating 
and evanescent components of both the input and the boundary Furthermore , 
both the propagat ing and evanescent components of the reflections from the 
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boundary are altered by the presence of the near field component of the input. 

This modified response of the boundary is included in the analysis as outlined 
by Appendix III. It is assumed, however, that the flow f ield produced by the line 
force input is not altered by the presence of the boundary . (In practice, this 
would not be the case since the impedance "seen" by the input would change due to 
the close proximity of the boundary. ) This simplification can be thought of as 
an alteration of the forcing function such that the characteristics of the flow 

field associated with the input remain the same 

Figure 17 shows that intensity estimates of the FREE FIELD method are nega- 
tive (and therefore unusable) below 400 Hz. Interestingly, the accuracy of the 
FREE FI ..D method actually improves in the 400-900 Hz range (compared to figure 
16) due to the close proximity of the source. Above 900 Hz, the FREE FIELD es- 
timates rapidly degenerate as they approach the secondary impedance mismatch. 

In contrast . the DIRECT method is seen to be accurate over the entire frequency 
range Thus , it may be concluded that the magnitude and character of the near 
field error associated with the FREE FIELD measurement method may also depend on 
the combined characteristics of the sources and boundaries in the vicinity of the 
measurement probe 

CONCLUDING REMARKS 

Overall , the results of the study indicate that the near field error asso- 
ciated with the FREE FIELD (two accelerometer) measurement method is much more 
serious than previously believed 

Simulations of FREE FIELD measurements near to simple sources show that 
the error curve of this method does not necessarily "hover" about the axis of 
0 error as implied by the more limited results presented in references 4, 5, 
and 6 Furthermore, it was shown that the frequency range of maximum error 
often coincides with the frequency range of maximum power flow when measuring 


„,.r simp!. source . I. - — «-* **» «» “ d th ° "* 5 " it “ d * "* 

yield error ...oci.t.d eith the FREE FIELD method I, very «~i«" *• *>“ 

distance b.t«..n th. measurement pr°»" and «h. sourc. And it ... .bo-n th.t th. 

accuracy of U. FREE FIELD method ia highly d.p.ndant on «h. radial and ■»-*•»> 

. . _ C vi t h respect to the source) when measuring 

locnt ion of the measurement probe (with re p 

near higher order enure., (dipole, end qu.drupol.s) . 

Simulation, of FREE FIELD m.a.ur.ment. ne.r to simple hound«i.e aho. that 

thi. method can actually indicate that the intensity is f lo.ing in the erong 
direction- Furthermore, it ... eho.n that this type .1 error i. not „.ce..arll, 
restricted to the loser I regency range, end can occur in frequency range, -her, 

the peer flo. i. relatively large and stable. 

Th . result . of tb. simulations lor tb. DIRECT method ere in .harp contrast to 

the results obtained lor tb. FREE FIELD method. In essentially all c.s.s und.r 

Investigation (.11 simple source, and boundaries studied, tb. direct Unit. 

difference approach produced an accurate estimate of tb. intensity fleeing in th. 

pl.t. Th. estimates produced by this method e.re so accurate . in fact . that no 

c found between the estimates and the exact solution for 
Bignif leant difference was found between 

any near field measurement condition studied. 

In fairness to the FREE FIELD approach it should be noted that many of the 
serious ne.r field errors inherent ,n the method can be suppressed for a given 

measurement condition by increasing th. spacing b.tv.en the tec acc.lenom.t.r. . 

Using this technique . th. near field .f « “* ^ 

frequency r.ng.s (beyond th. r.ng. of interest) This solution is not a pan... . 

since item also cause the error associated with the f init.dlff • enc. 
how«v«r, since n wm. 

approximations to greatly increase in significance 

u should also be noted that the FREE FIELD approach may be quite useful for 

f fhe nnwer f low m situations where there are no sources 
obtaining estimates of the power t 
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nearby , or in situations where the properties of the structure do not change in 
the vicinity of the measurement probe . 

Lastly , it should be pointed out that both methods are subject to other types 
of measurement error. Quinlin-* and Redman- White 15 have shown that the FREE FIELD 
method has difficulty producing an accurate estimate of the structural intensity 
in reverberant flow fields Furthermore , Hickol ' s results 6 suggest that the 
presence of the measurement probe itself can significantly alter the flow field 
and therefore degrade the accuracy of the measurements It is expected that the 
DIRECT method suffers from similar adverse effects . 


i 
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APPENDIX I 

EQUATIONS FOR THE TRANSVERSE VIBRATION OF AN INFINITE BEAM 


Genera] Solution for a Discrete Forcing Function 


The governing differential aquation for the transverse vibration of an 
infinite bean is given by (see reference 7) 




>1 > ... > 


1)1- 


(- 1 . 1 ! 


If it as suned that theforci ng f unct i on i s s imple harmonic in character and 
occurs at a discrete Locat ion in space, i . e . 


pl r , I • I r i, 


(.‘ 1.21 


then the result ing mot ion of the beam is also simple harmonic and the governing 
equation becomes 


m j. - I ' 

V 1 In r I, 

n n 


(•121) 


Now for the case of an infinite beam 
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. 1 1 


. ' r 


(. 1 . 1 ) 


and def ining I. 1 - II. the governing equat ion becomes 


-t I' 

I , I * 1-7 >1 r 1 1 ' • 

i,t I! 


II 


•l.M 


For values of r .• tl , the governing equation is 


I , 1 - M • I: - | ri 0 

■ Jr 1 


.1.11 1 
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This linear ordinary homogeneous differential equation is easily solved using 
operator methods and the solution is given by 

0 .'('if;*' t C>'J kr i ('r + (V" *"]< J+ ,( , I A. 

for values of ,r / I). Furthermore , from physical considerations , the value of >/ 
must remain finite as r approaches infinity. Therefore, C 4 0 and since all 
waves propagate away from the source, ( \ I) for values of > 0. Thus , the 
solution becomes 


for values of .r ■ f) 




Solut ion for a Point Force Input 


For the case of a discrete point force of amplitude located at ,r 0 on an 
lnfinitebeam, the following two boundary conditions apply : 


l.r (II 0, 


/((>, /’ 

l>i ' *»» H >.r II) 

. > / 1 ■> 


Applying the boundary condition of equation (A . 9) to the solution given b; 


equation (A 8 ) , it car. be shown that 


f i il ’i 


Thur the solut ion reduces t< 


i / 1 i . M f 




Applying the boundary condition of equation (A . 10) to the solution given by 


(A 12) , it is found that 


d W) 




and the solution for the transverse displacement for .r ■ I) as a function of time 
is given by 


*'•" umV' '*' " 


[.4.1 It 


Differentiating with respect to time to obtain the transverse velocity for 


<1 , 


nl.r.t ' | 


I \Iih ') 


>»T j , 


(.4.1.1) 


The angular velocity for j ■ I) is found by differential ing equation (A . 15 ) 
with respect to ./■ Thus 


;>r t l/tt- ) 


f L i f ii' 1 , 


(.4 in) 


The moment in the beam for r 0 as a function of time can be obtained by 
differentiating equation ( A 14 ) twicewith respect to j . The result i s 


Wil-r.n H 


J 


i - l !( 


^ r . 1 r ). I 


The shear in thebe am fere be an similarly be ob tamed f rom the third order 
differentiation of (A 14) with respect to < as follows: 


/■' 

t.n . a : • '*■ 

■I r ' I 


( .1.1 4) 



. , -V 


Solution for a Point Moment Input 


For the case of a discrete point moment of amplitude M,, located «t / - 0 on 
an infinite beam, the following two boundary conditions apply: 


i/(x - 0) • 0, M-W 

M h {r - 0 ) B^’hr-O) ~ Mp/2 . [A.2U 

ax* 

Applying the boundary condition of equat ion (A . 19) to the solution given by 
equation (A . 8) , it can be shown that 


Cl C, 


Thus the solution reduces to 


1 /l.r.M (',!-• > Kr (.'1.22 

Applying the boundary condition of equat ion ( A . 20) to the solution given by 
( A . 22 ) , it is found that 


1 M /#**)' 

and the solution for the transverse displacement for x It as a function of time 
is given by 


f/l x./ 1 


■ 'r i i j , kr , li." 

i in-) 


Differentiating with respect to time to obtain the transverse velocity for 


r II , 
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Th« *ng„i„ velocity lor r , 0 i. f by dif f 
with respect to r , Thus 


ng equation (A . 25) 


<>r (4 JU-) 




Th« moment in the beam for ,r 


— • U as a funct ion of time can be obtained by 

differentiating aquation (A 24) t.ie. „i,h ra.p.ct to r Th. r.ault i. 


-> Mr./) ttifv :K 

(Kr- I 


If J Kr a f 


The shear in the beam for , 0 can s imi larly be obtained f rom th# third order 

di f f wnti.t ior of (A 24) »i th respect to , as follows : 


1 1 h/ 

Q\*.n JJ i 

/hr 1 


,trr t f 



* 

L 


APPENDIX II 


EQUATIONS FOR THE TRANSVERSE VIBRATION OF AN INFINITE PLATE 


Governing Differential Equation for a Pvint Force 
The governing differential equation for the trenaverse vibration of an 
infinite plate is given by (see reference 7) 


it~n 

HV A i I t »i? , l> r,<t>J). 

(It ~ 


(A:w) 


If it assumed that the forcing function is a simple harmonic force tthich 
occurs perpendicular to tfcie plate at a discrete point in space , i . e . 


/'( v. O- 1 ) /’<M r h , 


/1.30» 


then the resulting motion of the plat e is also simple harmonic and the governing 
equation becomes 


H 






4.31 > 


Defining J'* mw - / / , this becomes 


r 


(V li 4 )v fi(r). - 


■1.321 


Fcr values of i ■ (I , the governing equat ion is 


(,1.33] 


Since a point force produces a symmetric response in the variable r , there can be 
no v dependence in the response i / Therefore , the Laplacian operator in this case 
reduces to 
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v2 <-> 2 . 1 d 

Hr* r Dr ' 


and tha biharmonic operator is given by 


V-V-'. M .:r, 

Substituting this lost result, the governing equation becomes 

(V- A- )( V J f k ’ ) >i (). (/t.M 

Solution for a Point Force 

Equation (A . 36) is a linear ordinary homogeneous differential equation. 
Using operator methods, this equation is easily solved. The solution is given by 


’! ■’) *'/'<> ;j ' r , 


where must satisfy 


,r- | ,i 
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and where >/' must satisfy 
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it) ’ r Hr 
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for va Lues of 


Equat ions (A 38) and (A 39) are alternate forms of the ordinary and modified 
Bessel’s equation of Oth order . respectively The solutions to these two equa- 
tions ars well known and are given , respect l vely , by 


>i f i i if i • < :, u f i, 
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V + jkr)-* (\H { 0 2 ) ( jkr). 


1 2 1 

where *(■ ) is the Oth order Hankel function of the 1st kind and Hq (■) is the Oth 
order Hankel function of the 2nd kind 

It can be shown (see ref erence 7) that Hq' *{ kr ) corresponds to an inward trav - 
oiling wave and that //q ' *( -jkr ) corresponds to a solution that grows exponen- 
tially without bound in kr . Therefore, from physical considerations, the coef- 
ficients of these two terms must be zero . Thus , the solution becomes 


\C 2 H { 0 2 '(kr) +■ CUH^i-jkr)^. 


(.4.42 i 


where H^ikr) corresponds to an outward travelling wave and i/g 2> ( jkr) corre- 
sponds to an exponently decaying solution in At. 

Applying the boundary condition 


(from symmetry considerations) to the solution given by (A .42) , it is found that 
( 'l < '_> . Thus , the solution is reduced to 


i/I r. t ) r : |//, , r'(AT) //J 2, ( jkr)}')*'. ( -141 

Now from reference 7 , the equation for the shear m the plate is given by the 
equat ion 


Q i 

< tr 


Applying the boundary condition 
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it is shown in reference 7 that the value of ('2 is given by 


r }Fp 
2 (HPk 2 ) ' 


(.4.47) 


Thus , the final solution for the displacement of the point driven infinite plate 
is given by 


»/<M) - (g ^,l //J 2 W> 2) i-jkT)\v»*. (.4.48) 

Equations for transverse velocity , angular velocity, bending moment , and 
shear can be obtained, respectively, by straightforward differentiation of 
equation (A . 48) using the following formulae : 
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Carrying out the operations indicated in equations (A 49) through ( A . 52) , the 
equat ions for transverse velocity , angular velocity , bending moment , and ehear 
for the case of an infinite plate driven by a point force are given, respectively, 
by: 



//o~ *( ,;A-r )!»- 

7U.’/ 
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Solution for a Point Moment 


The solution for an inf inite plate driven by a simple harmonic point moment 
can be obtained directly from the solution for a point force by simple differen- 
tiation with respect to the cartesian coordinate x (see reference 7) . Thus, the 
solution for the transverse displacement due to a point moment is given by 




and from the calculus it is fcnown that 


i) si tup d 

rostp . 

dr r dtp 


Since the monopole solution contains no <t> dependence, the dipole (point moment) 


solution is given by 





Performing the indicated operat ion , the result is 

i, • jU\ l \ )kr )>>-•', M-W) 

Note that this solution for the transverse displacement of the plate depends on 
both the radial d 1 stance r , and the azimuthal angle <p . Therefore , the power flow 
field due to a point moment input is not radially symmetric, and the flow field 
contains a component due to twist in addition to the ordinary shear and bending 
component s 

Equat ions for transverse velocity, angular velocity , rata of twist , bending 
moment , shear , and twisting moment can be obtained, respectively , by straight- 
forward differentiation of equation (A 58) using the following formulae 
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Carrying out the operations indicated in equations (A. 59) through (A. 64), the 
equation* for transverse velocity, angular velocity, rate of twiet , bending 
moment , shear , and twisting moment for the case of on infinite plate driven by a 
point moment tire given, respectively, by: 
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Solution for a Lateral Quadrupole 

The solution for an infinite plate driven by a simple harmonic point lateral 
quadrnpole can be obtained directly from the solution for a point moment by simple 
differentiation with respect to the cartesian coordinate y Thus , the solution 
fot the transverse displacement due to a lateral quadrapole is given by 


HyiHBk) 1 1 


i'-M, 


\kr) 1 ]H\ 2 \ - jk-r )]^^ 1 , 


( 4 . 71 ) 


and from the calculus it is known that 




Equations for .mi,.,,. ,alocit y . angular ,. looity , „ t . „ ( _ b „ djng 

nomant. ahaar. and t .i at mg moaumt can ba obtain.d, raapactival,. b, .tr.ight- 
forward dif f urantiut ion of aquation (a. 75) using th. loruiula. gi,.„ by , qu . tion8 

(. 59, through (. M) . Parking opurutions mdicut.d in tb.a, on 

equation ( A . 73) . the results are given by 




11 I^W) f -jbr)}r^'. 


k “’ Q >‘ 


jH \» { I H?\ - jkr 


ro.s2<p 
(HR) I* 


&<P I >\ i *> \ 

r j i ^ l (kr) ^ // 2 { ' J fir )j , 


(>1.74) 
, (4.75) 
(-4 76) 


U J*‘ 2 Vp 

’lb ‘ I'D.ldwiit) ( ’“‘I 

r ( - I . 


■Wj (&r) t //*,“'( - jkr )| 


;* r ,| * ^[//i’W) . «5"< >*D| 

,//| 2l ( VM! £'p!»5 ,, (*r| 4 «;»( y *a| I 


I ^ I I / 1 - * 
■ (* 


(4.77) 






/< 

•» '.s <,*>.<! i ij rp # J' 


"a[!//: 31 


(‘■’ I * jHj 11 ! j*r)| 


2 
(/•r 




v, , 


II I 


Ji'-Q, 


r,M -* ''“'•[(Jl-r) l "5 2l(frr> > // ! J, ( 

jk >') I 


(4.78) 


< 4.7ft) 


39 


APPENDIX III 


EQUATIONS FOR THE TRANSVERSE VIBRATION 
OF BEAMS NEAR TO SIMPLE BOUNDARIES 

Introduction 

This appendix is largely a recapitulation of the theory developed by B . R. 
Mace in reference 8 . For more detailed explanations of the concepts presented 
here see reference 8 

As shown in Appendix I , the solution for the transverse di sp lac ament in a 
beam due to a simple harmonic disturbance is given by equation (A. 7) (repeated 
here for the convenience of the reader) 


'/< '.*) <’v > kr • Cy '* J ■ (V kr t C 4 ,- k -r\ t 


l. J 


where the wavenumber k is given by k* m^-f H . It should be noted that the 
Coefficients in equation (A 7) may have complex values . Als ' note that < ] and < 
correspond to the propagating and decaying portions of a right traveling wave , 
respectively Simi larly , ( . and < ( correspond to the propagating and decaying 
portions of a left traveling wave, respectively. 

Now define the following vectors 



and also define the following matrices 



( A* 


r 



Utilizing equations ( A . 80 ) and (A . 81 ) , a vector describing the magnitude of the 
transverse displacement of the beam at any location x can be expressed as 


<H-r) ((/)(?* f f/j '<7- ), 


1,1.8 


and at the position .r 0 this vector is given by 


VlO.n i [7 j,7 ). 


( , 1.8 


where [ / 1 is the ident i ty matrix . 

Utilizing equation (A 82) , define the following vec cors : 


" '(■r) |/i<7 f ; .7 (.»■) i/l - 1 , i . 


( ,4.8‘ 


Using equation (A 84) it can be shown that the amplitudes of the displacement . 
slope , bending moment , and shear of the beam at any location x can be computed 
using the following matrix equation (see reference 8) : 
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Wave Ref lection at Point Discontinuities 

Mow assume that the location .r 0 coincides with the location of a point 
discontinuity (boundary) which is characterized by the following dimensionless 
translational and rotational stiffness parameters , respectively : 


A,, ( u- ’m, *- _/u.2^ t) ( *- A*,, )/( Elk * ) , ( .4 .8<> ) 

l\ it I hptno)'' ! k'o ) ( El k ) , (.4.87) 


Further suppose that a set of right traveling waves represented by u* ‘ is 
incident on this boundary and gives rise to a set of reflected (left traveling) 
uaves represented by •( given by the equation 


( , 4 . 88 ) 


where ' i is the reflection mat r ix . 
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It i. then in r.f.r.nc.8 th.t th. r.flhetiOhMtrU |r|. for a boundary 
ch.r.ct.rU.db, aquations U 86) and CA.87) , i, e iv.»b, th. aquation 


"’(I l) +P '(/ »)' 


wh®r« 


r ” T: u -"( t j ’ 


(A. wo I 


14 *- ( 1 j)hfl i 


Sub.Utut.n6.qu.tten (. 89, into .qu.t.cn <*.80 . th. v.ctot. tor th. 
r .ght tr.v.ling and l.ft tr.v.ling eo^pon.nt. of th. ..«»«■• di.pl.vo.nt .t 
location on th. .nc.d.nt s.d. of th. b.« ( r 0) are giv.n b. 


,1 ■ ( r I : f'-'i' ;*» <■ 


Tb. v actor .qu.t g.v.n b, (A .92) can no. b. .ub.t.tut.d into th. voter .qua- 
t , on of (A .85) Thu. . .qu.t.on (A 85) can b. u..d to c.lcul.t. th. tran.vrv 




m 


d ” Pl * C, "* nt ' ,h * "I"- *>” >«««•—«.- the ahn.ro, th. b..a, t lo . 
cation on th* incident side of tha baajn (r 0). Th. tr^a.ar,. „„ „ gul „ „ loe . 

*' 1 *’ °‘ th * b *“' 6 ' <5b, * i - d '«■ “>• tranavarsa diapl.c.^t an d , lop , 

of tha baam, raapactivaly, by di„.ranti.tio«.ithr.. p .et to ti«a. 
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